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Abstract 

A search for pair-produced leptoquarks has been performed using a sample of e + e~ colli- 
sion events collected by the OPAL detector at LEP at e + e~ centre-of-mass energies of about 
183 GeV. The data sample corresponds to an integrated luminosity of 55.9 pb _1 . The lepto- 
quarks were assumed to be produced via couplings to the photon and the Z° and then to decay 
within a single fermion generation. No evidence for contributions from leptoquark pair produc- 
tion processes was observed. Lower limits on scalar and vector leptoquark masses are obtained. 
The existing limits are improved in the region of large decay branching ratio to quark- neutrino. 
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1 Introduction 



In the Standard Model (SM) quarks and leptons appear as formally independent components. 
However, they show an apparent symmetry with respect to the family and multiplet structure of 
the electroweak interactions. It seems therefore natural that some theories beyond the SM 
predict the existence of new bosonic fields, called leptoquarks (LQs), mediating interactions 
between quarks and leptons. The interactions of LQs with the known particles are usually de- 
scribed by an effective Lagrangian that satisfies the requirement of baryon and lepton number 
conservation and respects the SU(3)c <S> SU(2)l ® U(l)y symmetry of the SM 0. This results 
in nine scalar (S) and nine vector (V) leptoquarks, grouped into weak isospin triplets (Si and 
Vi), doublets (5*1/2, §1/2, V1/2 and V1/2) and singlets (So, S , V and V ) []. They are shown in 
Tables [I] and |^. Under these assumptions, only the mass and the couplings to right-handed 
and/or left-handed fermions, denoted by Ar and Al, remain as free parameters, since the cou- 
plings to the electroweak gauge bosons are completely determined by the electric charge and 
the third component of the weak isospin. Moreover, in order not to contradict the existing 
indirect constraints on the leptoquark masses and couplings coming from low energy data such 
as rare decays [@,||[], the requirement that a given LQ couples to just one family of fermions is 
imposed. 

Several experimental results constrain the existence of leptoquarks. Searches for events 
with LQ single production, where a first generation LQ could be formed as a resonance be- 
tween an electron and a quark, were performed by the ZEUS and HI experiments at the ep 
collider HERA || and by the LEP experiments ||. In e + e~ collisions the quark comes from a 
resolved photon emitted by one of the LEP beams. As the production process directly involves 
a LQ-lepton-quark interaction, limits on the LQ mass, Mlq, can be derived as a function of 
the couplings, A, to fermions. Leptoquark masses below about 80 GeV are excluded for A 
values greater than a few 10~ 2 (about one order of magnitude smaller than the electromagnetic 
coupling A e — 0.3). For A = A e , HI excludes LQ masses up to 275 GeV. Both LEP and 
FERMILAB experiments have searched for events with LQ pair production 0,0], setting limits 
on Mlq as a function of ft, the branching ratio of decay into a charged lepton and a quark. 
On the contrary, these limits do not depend on the couplings to fermions because the LQ pair 
should be produced via coupling to the gauge bosons. First generation scalar LQs with masses 
lower than about 200 GeV are excluded, assuming ft > 0.5, while for (3 = the mass limit is 
79 GeV (DO). Second generation LQs are excluded below 160 GeV, assuming ft = 0.5, and 
below about 200 GeV, if ft = 1. Third generation scalar leptoquarks with charge |Q e . m .| = § 
or I and ft = 1 are excluded for masses lower than 99 GeV by CDF, while DO excludes third 
generation scalar leptoquarks with charge |Q e . m .| = | and ft = below masses of 94 GeV. 
In the same cases vector LQs belonging to the third generation are excluded for masses below 
170 GeV (CDF) and 148 GeV (DO). 

In principle, LQs of all three generations can be pair-produced in e + e~ collisions at LEP, 
by 5-channel 7 or Z° exchange and, in the case of first generation LQs, by exchange of a quark 
in the i-channel [|J. Due to the existing upper limits on the couplings, A, to fermions, the 
i-channel contribution to the first generation production cross-section is negligible so that, for 
a given LQ state 0, the cross-section depends on the mass, the electric charge and third com- 

: In this paper the notation used in j| is adopted. This is slightly different from the notation used in 0. 
2 In this paper "state" will be used to refer to a charge eigenstate within a multiplet. It will be denoted by 
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Table 1: Quantum numbers and couplings for scalar leptoquarks. Q e . m . is the electric charge 
in units of e, 1% the third component of the weak isospin and f3 denotes the branching ratio 
of decay to a charged lepton and a quark. Under the assumption of coupling within a single 
generation the same table must be repeated for second and third generation with the obvious 
substitutions e — > /i, r, u — > c,t and d — > s,b. 

ponent of the weak isospin, but is independent of the A couplings. On the other hand, for 
couplings smaller than O(10~ 5 ) the lifetime of leptoquarks would be sufficiently long to pro- 
duce a secondary decay vertex, clearly detatched from the primary vertex, or even outside the 
detector. This topology is not considered here as the charged tracks are required to come from 
the interaction vertex. To summarize, the present analysis covers the region of values of the 
couplings to fermions from about 1CT 6 to 1CT 2 . The decay of a heavy LQ into a quark and 

Si (Qe.m.) or V 1 (Qe.m.), where the index I represents the weak isospin while the number in parenthesis is the 
electric charge in units of e. 



5 



LQ 


h 


Qe.m. 


decay 


coupling \l,r 





V 





-2/3 


eld R 
e~ R d L 
v e u R 


^RV 
AlVo 


2X LV +X RV 


V 





-5/3 


e R u L 


'Wo 


1 


v 1 


1 



-1 


1/3 
-2/3 < 
-5/3 


vJ-r 

V e U_R 

[ e ldR 
eluR 


V 2A i y 1 
— AlVi 




1/2 

1 


v 1/2 


1/2 
-1/2 


-1/3 

-4/3 < 


V e dR 
. e RUL 

e L d R 

[ e R d L 


Alv 1/2 
Aw 1/2 

Alv 1/2 

V/2 


X 2 
RV l/2 

«l/2 + ^1/2 

1 


Vl/2 


1/2 
-1/2 


2/3 
-1/3 


V e UR 

elu R 


A iVi /2 




1 



Table 2: Quantum numbers and couplings for vector leptoquarks. Q e . m . is the electric charge 
in units of e, I3 the third component of the weak isospin and (3 denotes the branching ratio 
of decay to a charged lepton and a quark. Under the assumption of coupling within a single 
generation the same table must be repeated for second and third generation with the obvious 
substitutions e — > /i,t, u — > c,t and d — > s,b. 

a charged lepton leads to final states characterized by an isolated energetic charged lepton, 
while for decays into a quark and a neutrino, the final state will have large missing energy. 
Therefore the following topologies can be considered for events that result from the decay of a 
leptoquark-antileptoquark pair: 

Class A: two hadronic jets and two neutrinos; it includes the final states u e T7^uu, VeT^dd, u^TJ^cc, 
v,J7^ss and v T V^\jb. 

Class B: two hadronic jets, one neutrino and one charged lepton of the first (u e eud) or second 
generation (u^/ics). In the hypothesis that each LQ couples to just one family of fermions, the 
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topology including two jets, one neutrino and one r lepton is not possible at LEP, since one of 
the two jets would have originated in a top-quark. 

Class C: two hadronic jets and one pair of oppositely charged leptons of the first or second 
generation, for example e + e~uu and /j, + fi~cc respectively. 

Class D: two hadronic jets and one pair of oppositely charged r leptons, r + r~bb. This case is 
considered separately from Class C because of the different possible r decays. 

In this paper a search is presented for pair-produced scalar leptoquarks of all three genera- 
tions performed with the OPAL detector. Additionally, at this time, it was extended to vector 
leptoquarks of the first and second generation. The pair production process gives the advan- 
tage, with respect to single production by electron-quark interactions, that all states can be 
produced, included LQs coupling exclusively to neutrinos. Therefore searches for this channel 
at LEP give the possibility to explore the region of large decay branching ratio into quark- 
neutrino, where the FERMILAB experiments have reduced sensitivity. Moreover this is the 
first search for pair produced vector LQs at LEP. The study is based on data recorded during 
the 1997 LEP run at centre-of-mass energies, y/s, between 181 and 184 GeV, corresponding to 
an integrated luminosity of 55.9 pb _1 . The luminosity weighted average of yfs is 182.7 GeV. 
Limits are derived under the assumption that only one state contributes to the cross-section. 
In the case of the third generation, mass limits will be given for scalar LQs which decay into a 
b-quark and either a r-lepton or a r-neutrino. 



2 The OPAL Detector 

The OPAL detector is described in detail in Ref. |5j. It is a multi-purpose apparatus having 
nearly complete solid angle coverage 0. The central detector consists of two layers of silicon 



micro-strip detectors [lOj surrounding the beam-pipe and a system of tracking chambers inside 
a 0.435 T solenoidal magnetic field. This system consists of a high-precision drift chamber, a 
large-volume jet chamber and a set of z-chambers measuring the track coordinates along the 
beam direction. A lead-glass electromagnetic calorimeter is located outside the magnet coil and 
covers the full azimuthal range in the polar angle range of 

| cos# | < 0.984. It is divided into two regions: the barrel (| cos# | < 0.82) and the endcaps 
(| cos 6 | > 0.81). The magnet return yoke, divided into barrel and endcap sections along with 
pole tips, is instrumented for hadron calorimetry in the region | cos# | < 0.99. Four layers 
of muon chambers cover the outside of the hadron calorimeter. Close to the beam axis the 



forward calorimeter and gamma catcher together with the silicon-tungsten luminometer [ITI 
complete the geometrical acceptance down to 33 mrad from the beam direction. 



3 The right-handed coordinate system is defined so that z is the coordinate parallel to the e + and e~ beams, 
with positive direction along the e~ beam; r is the coordinate normal to the beam axis, <p is the azimuthal angle 
with respect to the positive direction of the x-axis (pointing towards the centre of LEP) and 9 is the polar angle 
with respect to +z. 
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3 Monte Carlo simulations 



Neglecting the t-channel contribution, the differential cross sections for the production of a pair 
of leptoquarks of mass Mlq in e + e~ collisions at a centre-of-mass energy y/s are || : 

W f l ~ 4M LQ / S ) f sin 2 9 \Us)\ 2 



das 
dcos 6 



8-s- 



day 
dcos 9 



3vra 2 



(l - 4M 2 LQ /s 



3 1 



1 + 



a=L,R 

4Mj> Q / S 



sm 



a=L,R 



for scalar and vector LQs respectively, a is the electromagnetic coupling and 



(1) 
(2) 



ka(s) 



Qe.m. Qa (^) 



S — M§ + iM Z T Z 



Q Z (LQ) 



(3) 



where Q e . m . is the LQ electric charge, M z and T z are the mass and the width of the weak 
neutral current gauge boson, and the couplings are given by 

Q Z (LQ) = 13 " ^.m.^6 W 



Qf(e) 



cos sin 9\y 

-\ + s\n 2 e w 

cos sin 9\y 
tan 



(4) 



In equations f| J3 is the third component of the LQ weak isospin and 9w is the Weinberg angle. 



The Monte Carlo generator LQ2 [|12| was used to simulate leptoquark pair events. Initial 
state QED radiation was included. The leptoquarks are assumed to decay isotropically in their 



rest frame and the hadronization of the final state qq' pair was performed by JETSET [13 



Samples of 1000 signal events corresponding to different values of the leptoquark mass, from 
Mlq = 50 GeV to M L q = 90 GeV in steps of 5 GeV, were generated for all the different 



decay topologies. The full simulation of the response of the OPAL detector |14| was performed 
on the generated events. Since they carry colour, leptoquarks may hadronize before decaying, 
if the couplings to fermions are small. This possibility was taken into account by introducing a 
systematic error on the detection efficiencies which was evaluated by using Monte Carlo samples 
in which the leptoquarks hadronized before decaying. 



All relevant Standard Model background processes were studied using Monte Carlo gener- 
ators. Two-fermion events (Z°*/7* — > fft/y), f = q,r) were simulated with PYTHIA [13| and 
KORALZ |g . The Monte Carlo programs HERWIG |T§ and PHOJET were used to 
generate two-photon hadronic events. Other processes with four fermions in the final state, 



including W pair production, were simulated with grc4f |T8[ and Vermaseren |T9] 



4 Analysis 

Charged tracks used in the calculation of physical variables were required to have their origin 
close to the e + e~ interaction point, to have at least 20 measured space points in the jet cham- 
ber and at least 50% of the hits geometrically expected. The minimum transverse momentum 
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of the track with respect to the beam direction had to be greater than 120 MeV. Electro- 
magnetic clusters were required to have an energy of at least 100 MeV in the barrel and 250 
MeV in the endcaps. Endcap clusters were also required to contain at least two adjacent lead 
glass blocks. Clusters in the hadron calorimeter were rejected if their energy was smaller than 
0.6 GeV (2 GeV in the hadron poletips, i.e. for | cos# | > 0.91). To avoid double counting, 
calculations of experimental quantities such as visible energy, transverse momentum, etc., were 
performed following the method explained in [ED . 



All the different topologies of signal events (classes A to D) are characterized by large 
charged multiplicities and large energy depositions due to the hadronization of the qq' pair. 
In this analysis no attempt was made to identify the flavour of the quarks. Electron and 
muon identification, required in the selection of events of classes B and C, was performed 



by making a logical "OR" of different standard algorithms pl| . The electron identification is 
based on the energy match between a track and the associated cluster in the electromagnetic 
calorimeter and uses a minimum of subdetectors in order to give high efficiency, while the muon 
identification requires a minimum number of associated hits in either the muon chambers or 
hadron calorimeter strips. The energy of identified electrons was taken from the energy of the 
electromagnetic calorimeter cluster associated to the identified electron track. 



4.1 The jet jet v v channel (class A) 

Signal events of class A are characterized by an acoplanar pair of hadronic jets and missing 
energy. A few preselection requirements were applied to the data. To reduce beam-wall and 
beam-gas interactions, the fraction of charged tracks that satisfied the quality criteria given 
above was required to be greater than 0.2. A similar requirement was applied to the non- 
associated electromagnetic clusters. Both the number of accepted charged tracks and the 
number of accepted non-associated electromagnetic clusters had to exceed four. Finally, the 
total visible energy, E vis , was required to be greater than 0.25^ and smaller than 1.25\/s. 
After the preselection, the following cuts were applied to the data: 

Al) The events contained no identified charged electron or muon of energy greater than 
0.12^5. 

A2) A two-dimensional cut was made in the plane (Pt mss /v^) vs - I cos ($miss) |, where p™ lss is 
the missing transverse momentum with respect to the z-axis and cos(# m i SS ) the cosine of 
the angle between the missing momentum and the z-axis. The cut is shown in Figure [I]. 

A3) The events were reconstructed into two jets using the Durham algorithm f22fl : the sum 
of the energies of the two jets, Ej e t s , had to be such that 0.25a/s < Ej e t s < 0.75a/s. 

A4) The jets were required to be acolinear by asking that cos(#jj) > -0.1, where 9^ is the angle 
between the two jet directions. 

Cuts Al and A2 are useful, in particular, in rejecting two-photon events and greatly reducing 
Z°*/7* _> ff(j) background. Cuts A3 and A4 completely reject Z°*/7* —> ff(7) events, and 
reduce four-fermion background. 
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Data 


Background 


4-lermions 


77 


Z fj -> 11(7) 




{ro) 


(Al) 


16281 


14913 


685 


9308 


4920 


90.0 


88.4 


(A2) 


87 


80.9 


52.2 


1.6 


27.1 


59.4 


58.7 


(A3) 


53 


46.8 


42.1 


0.9 


3.8 


58.5 


56.9 


(A4) 


4 


3.6 


2.7 


0.8 


<0.1 


41.7 


39.8 



Table 3: The remaining numbers of events after each cut of selection A for various background 
processes are compared with the data. The background is normalised to an integrated lumi- 
nosity of 55.9 pb -1 . The last two columns report the signal efficiencies, in percent, for events 
with M L q = 85 GeV, for scalar and vector LQs respectively. 



Table |^ shows the number of events after each cut, compared with the number of background 
events as predicted from Monte Carlo samples, and the efficiencies for signal events correspond- 
ing to Mlq = 85 GeV. Four events in the data survive the selections, in good agreement with 
the number of expected background events from Standard Model processes, 3.57 ± 0.83 (stat.). 

4.2 The jet jet I ^ v channel (class B) 

After the application of the same preselection cuts described in section 4.1, the selection of 
signal events of class B proceeded as follows: 

Bl) The event was required to contain at least one identified electron or muon with an energy 
of at least 0.1 \fs. 

B2) No charged tracks and at most one (no) electromagnetic cluster within a cone of half- 
aperture 15° around the most energetic electron (muon). 

B3) The missing transverse momentum had to be greater than 0.09^. 

B4) A 2-dimensional cut in the plane M^ vs. cos(#i„) was applied as shown in Figure |2], where 
9\ u is the angle between the momentum vector of the most energetic lepton and the missing 
momentum vector, while Mh, is their invariant mass. 

B5) After having removed the most energetic lepton, the event was forced into two jets using 
the Durham algorithm. The sum of the energies of the two jets, Ej e ts, had to be such that 
0.25a/s < Ej Cts < 0.75>/i. Moreover, a 2-dimensional cut in the plane cos(#jj) vs. Mjj was 
applied as shown in Figure |^, where 6^ is the angle between the two jet directions and 
Mjj is the invariant mass of the two-jet system. 

Two-photon events contribute negligibly to the background after selection B3. Cuts B2 and 
B3 are particularly efficient in reducing Z°*/7* — > ff(7) events. 

In Table |] the numbers of events after each cut are shown, compared with the numbers of pre- 
dicted background events and the efficiencies for signal events corresponding to Mlq = 85 GeV. 
One event in the data is retained after the selection for the first generation and three events for 
the second generation, in good agreement with an expected background of 2.22 ± 0.21 (stat.) 
and 2.84 ± 0.18 (stat.) events, respectively. 
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Data 


Background 


A £ 

4-iermions 


77 


Z u /7 -> ff(7) 


B / O / \ 

e£ (%) 


B ( ()/ \ 

£v (%) 


(Bl) 


1110/304 


1106/321 


277/151 


28.4/2.5 


801/169 


89.1/87.7 


90.7/90.1 


(B2) 


196/109 


179/101 


128/99.6 


20.8/0. 


30.2/1.4 


79.8/80.3 


81.7/79.9 


(B3) 


108/95 


98.4/85.9 


94.6/85.5 


0.8/0. 


3.0/0.4 


73.5/72.3 


75.2/73.6 


(B4) 


23/29 


19.4/20.7 


17.1/20.3 


0.8/0. 


1.5/0.4 


59.8/58.8 


60.4/61.7 


(B5) 


1/3 


2.2/2.9 


2.0/2.6 


o./o. 


0.2/0.3 


43.1/49.1 


43.6/51.5 



Table 4: Tie remaining numbers of events after each cut of selection B for various background 
processes are compared with the data. The background is normalised to an integrated luminos- 
ity of 55.9 pb _1 . The last two columns report the signal efficiencies, in percent, for events with 
Mlq = 85 GeV, for scalar and vector LQs respectively. When two numbers are separated by a 
slash, the first one refers to the first generation and the second one to the second generation. 



4.3 The jet jet l + I channel (class C) 

Signal events of this type have small missing energy and are characterized by the presence of 
a pair of high energy charged leptons of the same generation that tend to be isolated. The 
same preselections as for classes A and B were applied. The following cuts were then applied 
to select events: 

CI) The presence of at least one pair of oppositely-charged identified electrons or muons was 
required. The most energetic leptons of the same generation and opposite in charge will 
be called the "pair" in the following. The energy of the most energetic lepton of the 
pair had to exceed 0.15\/s, while an energy of at least 0.1 y/s was required for the other 
lepton. For the first generation, the distributions of the energy of the most energetic 
electron of the pair for data, expected background and signal events corresponding to 
Mlq = 85 GeV, are shown in Figure |. 

C2) No charged tracks and at most one (no) electromagnetic cluster within a cone of half- 
aperture 15° around the most energetic electron (muon). The same requirements, but 
considering a cone of half-aperture 10°, were applied for the less energetic lepton of the 
pair. 

C3) The two leptons of the pair were required to be acolinear: cos(#u) > -0.85, where 6\\ is 
the angle between the two leptons. 





Data 


Background 


4-fermions 


77 


Z0*/ 7 * _> ff( 7 ) 


e°s (%) 


4 (%) 


(CI) 


97/10 


96.0/7.4 


19.9/4.3 


0.4/0. 


75.7/3.1 


80.7/80.3 


79.0/77.4 


(C2) 


2/2 


2.8/1.5 


2.7/1.5 


0.1/0. 


<0.1/0. 


67.2/69.2 


67.4/67.3 


(C3) 


2/2 


1.6/1.0 


1.6/1.0 


o./o. 


<0.1/0. 


63.7/65.7 


63.7/63.9 



Table 5: The remaining numbers of events after each cut of selection C for various background 
processes are compared with the data. The background is normalised to an integrated luminos- 
ity of 55.9 pb _1 . The last two columns report the signal efficiencies, in percent, for events with 
Mlq = 85 GeV, for scalar and vector LQs respectively. When two numbers are separated by a 
slash, the first one refers to the first generation and the second one to the second generation. 



Background from two-photon events is removed after selection CI. Cut C2 suppresses Z°*/7* — > ff(7) 
events. The lepton acolinearity requirement C3 is useful in further reducing four-fermion back- 
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ground. 



The numbers of events after each cut, compared with the numbers of expected background 
events and the efficiencies for signal events corresponding to Mlq = 85 GeV, are shown in 
Table |5|. After all cuts two candidates survive the selection both for the first and the second 
generation; 1.63 ± 0.14 (stat.) background events for the first generation and 0.98 ± 0.15 (stat.) 
for the second generation are expected. 



4.4 The jet jet r + r channel (class D) 

Signal events of this type are characterized by a pair of isolated r-leptons and a pair of energetic 
hadronic jets. The background comes predominantly from Z°*/7* — > ff(7) and four-fermion 
processes. The selection begins with the identification of r lepton candidates, which is identical 



to that in P3[ , using three specific algorithms to identify semileptonic and hadronic r-lepton 
decays. On average, 2.3 r candidates per signal event are thus identified. The original r-lepton 
direction is approximated by that of the visible decay products. The following requirements 
are then imposed: 

Dl) Events are required to contain at least nine charged tracks, and must have at least two 
r-lepton candidates, including at least one pair whose members each have electric charge 
|g| = 1 and whose charges sum to zero. Pairs not fulfilling these requirements are not 
considered further. 

D2) Events must have no more than 20 GeV of energy deposited in the forward calorimeter, 
gamma- catcher, and silicon-tungsten luminometer; a missing momentum vector satisfying 
| cos# m i SS | < 0.97, a total vector transverse momentum of at least 0.02 y/s, and a scalar 
sum of all track and cluster transverse momenta larger than 40 GeV. 

D3) Events must contain at least three jets, including single electrons and muons from r-lepton 
decay which are allowed to be recognised as low-multiplicity jets, reconstructed using the 
cone algorithm as in [2^] and no energetic isolated photons. An energetic isolated photon 



is defined as an electromagnetic cluster with energy larger than 15 GeV and no track 
within a cone of 30° half-angle. 

D4) Events must contain no track or cluster with energy exceeding 0.3-^/s. 

For events surviving these requirements, the tracks and clusters not belonging to the r pair 
(henceforth referred to as the "rest of the event" or RoE), are then split into two jets using the 
Durham algorithm. Two pairing schemes between the two r candidates and the jets are thus 
possible. The invariant masses m T j of the two resulting r-jet systems within each pairing scheme 
are then calculated using only the r-lepton and jet momentum directions and requiring energy 
and momentum conservation. The pairing scheme exhibiting the lesser difference between m T ji 
and m T j2 is then chosen. Then, in order for a r candidate pair to be considered further, the 
following requirements on m T ji and m T j2 are imposed, consistent with the hypothesis of the 
decay of two heavy objects of identical mass: 

D5) Both m T ji and m T j2 must be at least 30 GeV. 

D6) The difference in invariant masses must be no more than 30% of their sum. 
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The distribution of \m T ji— m T j2\/\m T ji+m T j 2 \ is shown in Fig. | (a) for the data, the background 
events, and for a signal sample with Mlq = 85 GeV. The resolution on m T j is typically less 
than 7 GeV, except very close to the kinematic limit. 



A likelihood method p5|,p6| is then applied to those events satisfying the above require- 
ments, in order to select a final r candidate pair for each event from those surviving, and to 
suppress further the remaining background. In each such event, for each remaining r candi- 
date pair and its associated hadronic RoE, a joint discriminating variable, £, is constructed 
using normalised reference distributions generated from Monte Carlo samples of signal and 
background events. The set of variables for the reference distributions includes some which 
characterize each of the two r-lepton candidates individually, some which describe their be- 
havior as a pair and some which characterize the RoE. For those variables describing the r 
candidates individually, separate reference distributions are generated for leptonic (electron or 
muon), hadronic 1-prong and hadronic 3-prong r candidates, in order to exploit the differences 
between the three categories. Distributions of some of these input variables as well as that of 
£ are shown in Fig. |5] (b) to (d). The r candidate pair having the highest value of £ is chosen 
as the definitive pair for each event, and the following requirement is then made: 

D7) £ > 0.93 

Table |] shows the numbers of observed and expected events after each requirement, along 
with the detection efficiency for a signal with Mlq = 85 GeV. Two events survive the se- 
lection while the background, predominantly from four-fermion processes, is estimated to be 
2.07 ± 0.15 (stat.) events. 





Data 


Background 


4-fermions 


77 


Z°*/ 7 * _> ff( 7 ) 


(%) 


(Dl) 


1322 


1070 


150 


791 


129 


59.4 


(D2) 


209 


191 


129 


3.0 


59.1 


57.8 


(D3) 


198 


181 


127 


2.9 


51.1 


57.5 


(D4) 


149 


139 


100 


2.2 


36.3 


56.2 


(D5) 


51 


55.5 


47.5 


0. 


8.0 


52.7 


(D6) 


41 


44.8 


38.1 


0. 


6.7 


50.2 


(D7) 


2 


2.1 


2.1 


0. 


<0.02 


32.6 



Table 6: The remaining numbers of events after each cut of selection D for various background 
processes are compared with the data. The background is normalised to an integrated lumi- 
nosity of 55.9 pb . The last column reports the signal efficiency, in percent, for events with 
M LQ = 85 GeV for scalar LQs. 



5 Results 

The detection efficiencies for the different topologies of signal events, as a function of the 
leptoquark mass M L q, are listed in Table [7]. 

The systematic uncertainties in the number of signal events for the following sources were 
evaluated (they are quoted as relative %): 

• The errors due to signal Monte Carlo statistics and the interpolation errors at an arbitrary 
point of Mlq were estimated to be 5-10%, depending on the signal topology. 
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M LQ (GeV) 


50 


60 


70 


80 


85 


90 


Signal topology 


Generation 














Class A (scalar) 


1,2,3 


14.3 


23.9 


31.4 


37.0 


41.7 


42.8 


Class A (vector) 


1,2 


15.0 


22.7 


32.3 


36.3 


39.8 


42.3 


Class B (scalar) 


1 


14.1 


25.5 


30.9 


40.0 


42.4 


45.2 


Class B (scalar) 


2 


15.1 


28.9 


38.2 


44.8 


49.1 


52.5 


Class B (vector) 


1 


13.9 


23.0 


35.5 


41.7 


43.6 


46.7 


Class B (vector) 


2 


15.3 


25.8 


38.5 


47.0 


51.5 


53.9 


Class C (scalar) 


1 


36.2 


45.5 


56.6 


59.3 


63.7 


64.5 


Class C (scalar) 


2 


37.1 


46.1 


55.9 


62.9 


65.7 


67.0 


Class C (vector) 


1 


33.8 


43.0 


52.8 


62.1 


63.7 


67.0 


Class C (vector) 


2 


35.5 


41.7 


53.4 


63.0 


63.9 


66.8 


Class D (scalar) 


3 


22.3 


27.4 


29.6 


32.0 


32.6 


32.8 



Table 7: The detection efficiencies for the various selections, in percent, as a function of the 
leptoquark mass, M L q. 



The error associated with the electron and muon identification method was evaluated 
using mixed events constructed by overlaying Z°*/7* — > qq events with single hemispheres 
of Z°*/7* — > 1 + 1~ (1 = e,/i) events at yfs ss M z o [P7| . Such events are topologically and 



kinematically analogous to qqlu events at y/s ~ 183 GeV. This error was found to be 
3.2% for electrons and 2.5% for muons. The error for r identification was evaluated to be 



1.2% [26 



The uncertainty introduced by the modelling of the variables used in the selections was 
estimated to contribute with 5-15%. This was evaluated by displacing the cut values 
by an amount corresponding to the difference between the mean values of the data and 
background distributions. The corresponding displacement in the detection efficiencies 
was taken as a systematic error. The contributions from each single cut were added in 
quadrature. 

The uncertainty in the flavour of the final state quarks contributes with less than 6% 
for class A and with less than 3% for class C. This was evaluated by comparing the 
efficiencies corresponding to all the different possible flavours in a given decay channel, 
characterized by the leptons in the final state (for example e + e~uu and e + e~dd for class C, 
first generation). The value of the efficiency was taken to be the mean value and the largest 
difference between the mean and the single contributions was taken as a systematic error. 
In the case of classes B and D the flavours of the two quarks are precisely determined 
by the hypothesis of coupling within a single generation and (for class D) the top-quark 
mass threshold and thus no error is assigned. 

The systematic error associated with the fragmentation of the leptoquark decay products 
was estimated to be about 10%. This was evaluated by comparing the signal efficiencies 
with the efficiencies obtained by applying the same selection cuts to signal events in which 
the two LQs fragment before they decay. 

The error on the integrated luminosity was 0.5%. 
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The errors were considered to be independent and added in quadrature. A total systematic 
error of 10-30% was estimated on the number of expected background events, from the fol- 
lowing sources: Monte Carlo statistics (10-20%), modelling of cut variables (10-25%), lepton 
identification (1.2-3.2%) and integrated luminosity (0.5%). 



No evidence of leptoquark production was observed in the data. Upper limits at 95% con- 
fidence level (C.L.) on the LQ pair production cross-section cx(e + e~ — > LQLQ) were computed 
from the observed numbers of events, the signal detection efficiencies and the number of ex- 
pected background events, using the procedure described in pg[] . The uncertainties both on 
the efficiencies and on the background were incorporated in the upper limits by numerical in- 
tegration, assuming Gaussian distributions, as suggested in |H| . To obtain the limit at a given 
point in the plane (3 vs Mlq, the independent analyses, corresponding to different classes of 
events, were combined by considering the total number of expected signal events. For example, 
for the first and the second generation, this number is given by 



N. 



tot 
exp 



a 



c 



2 e c + 2/3(1 - (3)e B + 



1-/3) 5 



(5) 



where a is the production cross-section, / £ is the integrated luminosity of the data and e A ' B ' c 
are the signal detection efficiencies for classes A, B and C respectively. 



Regions of the plane /3 vs Mlq excluded at the 95% C.L. were determined by comparing 
the upper limits on the production cross-sections with the total cross-sections computed by 
integration of Equations [l| and In the case of the third generation, limits have been evaluated 
for scalar states which decay only either into a r-lepton and a b-quark, or into a v T and a b- 
quark. For the third generation state S\i%{— 2/3) only (3=1 can be considered since below 
the top-quark threshold only the right-handed coupling is possible. Figures |6| (a),(b) and |7| (a) 
show the upper limits of the production cross-sections as functions of the LQ mass, for scalar 
LQs with (3 = and (3 = 1. The limit on Si/ 2 (l/3) is valid for all three generations but the 
same is not true for the limit on Si (2/3), since the third generation state, under the hypothesis 
of diagonal coupling, would decay to a top-quark and a v T . Figures || (a),(b) and || (a) show 
the upper limits of the production cross-sections as functions of the LQ mass, for vector LQs 
with (3 = 1, (3 = 0.5 and (3 = respectively. The regions excluded in the plane (3 vs Mlq of the 
states Si/ 2 (— 2/3), V\/2{— 1/3) and Vq(— 2/3), whose (3 can range from to 1 (0.5 to 1 for V ) 
according to the relative weights of the left and right A couplings, are reported in Figures ^ (b) 
and [5] (b) for the first and the second generation. The mass limits obtained from the present 
analysis are summarized in Table || Because of their much smaller cross-sections this search is 
not sensitive to the production of the states Sq(— 1/3) and Si (—1/3), so previous OPAL limits 
are quoted. 



6 Conclusions 

A data sample collected with the OPAL detector at y/s = 183 GeV, corresponding to an 
integrated luminosity of 55.9 pb _1 , was analysed to search for events with pair production of 
scalar and vector leptoquarks. The search included scalar states of all three generations and 
vector states of the first and second generation only. In the case of vector leptoquarks this is 
the first pair production search at LEP. The present analysis covers the region of small values 
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T ( \ 




P 


1 gen. 


2 gen. 


3 gen. 


Sq 


-1/3 


[0.5,1] 


44.2(*) 


44.2(*) 


— 


So 


-4/3 


1 


85.8 


85.5 


82.7 




-2/3 


[o,i] 


80.8(**) 


78.8(**) 


82. 2(* * *) 


Sl/2 














-5/3 


1 


87.0 


86.8 






+ 1/3 





71.6 


71.6 


71.6 


~ 

Sl/2 














-2/3 


1 


81.8 


81.5 


76.9 




+2/3 





84.8 


84.8 


— 


Si 


1 /o 

-1/3 


0.5 


A A O / ■ \ 

44. 2{*) 


A A O / . \ 

44.2(*J 






/I /Q 

-4/3 


1 


/ .0 


1 .0 




V 


-2/3 


[0.5,1] 


85.8(**) 


85.1(**) 


— 


V 


-5/3 


1 


90.5 


90.4 


— 




-1/3 


[0,1] 


88.0(**) 


87.5(**) 


— 


v 1/2 














-4/3 

/ 


1 


90.1 


90.0 






+2/3 





87.5 


87.5 


- 


V, In 














-1/3 


1 


88.8 


88.6 


- 




+ 1/3 





89.8 


89.8 






-2/3 


0.5 


85.8 


85.1 






-5/3 


1 


90.8 


90.7 





Table 8: Lower mass limits, in GeV, for scalar and vector leptoquarks as obtained from the 

present analysis. 

(* ) LEP1 limits from OPAL. 

(**) Minimum allowed value for Mlq. See Figures 7(b) and 9(b) for limits as functions of f3. 
(* * *J This limit is valid for [3 = 1 (i.e. Xl — below the top-quark threshold). 

of the couplings, A, to fermions (about from 10~ 6 to 10~ 2 ). No significant excess with respect to 
Standard Model predictions was found in the data. Lower mass limits for leptoquarks were set, 
under the assumption that only one leptoquark contributes to the cross-section. The results 
improve previous LEP lower limits on scalar leptoquark masses by 25-40 GeV, depending on 
the leptoquark quantum numbers. With respect to the existing limits from CDF and DO 
experiments, mass lower limits are improved by about 10 GeV, for first and second generation, 
in the region of small values of the branching ratio of decay to a charged lepton and a quark. 
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the jet jet v v channel: The distribution of events in the plane (p™ lss /\/s) vs. 
, shown after the preselection and cut Al, for the data (a), the simulated back- 
ground (b) and simulated scalar LQ signal events of class A with Mlq = 85 GeV (c). The 
area of each box is proportional to the logarithm of the number of events falling within a 
two dimensional bin and is normalized with respect to the total content of each histogram 
separately. 
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Figure 2: 

Class B, the jet jet I v channel: The distribution of events in the plane M Mi , vs. cos(0 MI/ ), 
shown after the preselection and cuts Bl - B3, for the data (a), the simulated background (b) 
and simulated scalar LQ signal events of class B, second generation, with Mlq = 85 GeV (c). 
The area of each box is proportional to the logarithm of the number of events falling within 
a two dimensional bin and is normalized with respect to the total content of each histogram 
separately. 
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Figure 3: 

Class B, the jet jet I v channel: The distribution of events in the plane My vs. cos(#jj), displayed 
after the preselection and cuts Bl - B4, for the data (a), the total simulated background (b) 
and simulated scalar LQ signal events of class B with Mlq = 85 GeV (c). The area of each 
box is proportional to the logarithm of the number of events falling within a two dimensional 
bin and is normalized with respect to the total content of each histogram separately. 
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Figure 4: 

Class C, the jet jet l + l~ channel: The distribution of the energy of the most energetic electron 
for events containing an electron pair after the preselections, for the data (points), the simulated 
background (full line, normalised to the integrated luminosity of the data) and a simulated signal 
(dashed line, arbitrary normalisation) corresponding to M LQ = 85 GeV. The arrow denotes the 
position of cut CI. 
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Figure 5: 

Class D, the jet jet t + t~ channel: Kinematic distributions for the data (points), estimated 
Standard Model background (full line, normalised to the integrated luminosity of the data), 
and a simulated signal (dashed line, arbitrary normalisation) corresponding to Mlq = 85 GeV. 
(a) Distribution of the difference in invariant mass of the tau-jet systems scaled by their sum 
after cut (D2). 

Figures (b) and (c) demonstrate the difference in the distributions of the same likelihood input 
variable for two different categories of r candidate, after cut (D4): (b) The momentum of 
leptonic r candidates, (c) The momentum of 1-prong hadronic r candidates, (d) The likelihood 
distribution after cut (D6). 
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Figure 6: 

(a) The 95% C.L. upper limits of the production cross-section for scalar LQs of the first and 
second generation with (3=1 (full lines) as a function of the LQ mass, compared with the 
theoretical production cross-sections, (b) Same as (a) for third generation scalar LQs. 
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Figure 7: 

(a) The 95% C.L. upper limit of the production cross-section for scalar LQs with f3 = (full 
line) as a function of the LQ mass, compared with the theoretical production cross-sections. 

(b) The region of the plane (3 vs M LQ excluded at the 95% C.L. for the state Si/ 2 (— 2/3) of the 
first and second generation. 
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Figure 8: 

(a) The 95% C.L. upper limits of the production cross-section for vector LQs of the first and 
second generation with (5 — 1 (full lines) as a function of the LQ mass, compared with the 
theoretical production cross-sections, (b) Same as (a) for vector LQs of the first and second 
generation with (3 = 1/2. 
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Figure 9: 

(a) The 95% C.L. upper limit of the production cross-section for vector LQs with (3 = (full 
line) as a function of the LQ mass, compared with the theoretical production cross-sections. 

(b) The region of the plane (3 vs M L q excluded at the 95% C.L. for the states V (— 2/3) and 
Vi/2(— 1/3) of the first and second generation. 
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